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Abstract 

The self-energy of the Dirac Brueckner-Hartree-Fock calculation in nuclear mat- 
ter is parametrized by introducing density-dependent coupling constants of isoscalar 
mesons in the relativistic Hartree-Fock (RHF) approach where isoscalar meson a , to 
and isovector meson it , p contributions are included. The RHF calculations with den- 
sity dependent coupling constants obtained in this way not only reproduce the nuclear 
matter saturation properties , but also provide the self-energy with an appropriate den- 
sity dependence. The relativistic density dependent Hartree-Fock (RDHF) approach 
contains the features of the relativistic G matrix and in the meantime simplifies the 
calculation. The ground state properties of spherical nuclei calculated in the RDHF are 
in good agreement with the experimental data. The contribution of isovector mesons 
TT and p , especially the contribution of the tensor coupling of p meson , are discussed 
in this paper. 



1 



1 Introduction 

It is well known that the relativistic mean field (RMF) theory has been extensively 
used to investigate the ground state properties of spherical and deformed nuclei [0, |^ in 
the past several years. Recently , it is used to study the properties of nuclei far from /3 
stability line^, ^. Bouyssy et al]^ extended the RMF to the relativistic Hartree-Fock(RHF) 
approach. The contributions of the exchange terms and isovector mesons in RHF to the 
ground state properties of nuclei have been emphasized. Since then, the RHF has been 
developed by several authors]^, 0]. 

Though the RMF or RHF approach has been quite successful in reproducing the bulk 
properties of nuclei. However, too large compressibility of nuclear matter is produced in these 
approaches . It might indicate an incorrect density and momentum behavior of the effective 
interaction described in the RMF and RHF approaches. Some discrepancies in reproducing 
the properties of finite nuclei in the standard RMF or RHF calculations are observed. The 
calculations yield correct binding energies but too little charge radius or vice versa. Such 
calculations reveal a new "Coester" band in a dependence of l/vch on Eb which is similar 
as in the nonrelativistic Brueckner-Hartree-Fock approach (BHF)[^]. A more fundamental 
and sophisticated way to deal with the many-body problem is the Dirac Brueckner-Hartree- 
Fock(DBHF) approach [9- 12]. Starting from a bare nucleon-nucleon interaction of one boson 
exchange potential, one solves the Brueckner-Goldstone equation in the nuclear medium. 
The Brueckner G matrix is the two nucleon effective interaction including nucleon-nucleon 
short range correlation. The DBHF has been quite successful in reproducing the nuclear 
matter saturation properties, as well as the compressibility |TT|, |T3|. It provides a correct 
description of the density- and momentum- dependence of the nucleon self-energy in the 
nuclear medium. However, the DBHF due to its complexity is mainly restricted to nuclear 
matter and only very few finite-nuclei calculations are performed so far[l^. Therefore, it 
is requested to remedy the deficiencies of the RMF or RHF approach with an effective 
interaction without losing the features of the relativistic G matrix and at same time retain 
the simplicity. 

Recently, there is a growing effort on developing the relativistic effective interaction both 
in nuclear matter and finite nuclei. Attempts have been also made to improve both binding 
energy and rms radii of finite nuclei simultaneously in the RMF or RHF with various effec- 
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tive interactions. Gmuca |]T5| has parametrized the results of the DBHF in nuclear matter 
in terms of the RMF with nonlinear scalar and vector self-interactions. Though the cou- 
pling constants obtained in such way is density-independent , these mesonic self-interactions 
implicitly represent the density- dependence. Brockmann and Toki |]16[ developed a relativis- 
tic density dependent Hartree (RDH) approach for finite nuclei. The coupling constants 
of isoscalar meson a and uj in the RMF are adjusted at each density by reproducing the 
nucleon self-energies in nuclear matter resulting from the DBHF instead of fitting them to 
the empirical nuclear matter saturation properties. The binding energies and rms radii of 
^^O and ^°Ca calculated with the density dependent interaction are in good agreement with 
experiments. However, Fritz, Miither and Machleidt [|T^] pointed out that the Fock terms 
are not negligible. The relativistic density dependent Hartree-Fork(RDHF) calculations were 
performed Due to the uncertainty of the DBHF at very low densities , the extrapo- 
lation procedure has to be adopted when applied to the calculation of finite nuclei. The 
sensitivity of the results calculated by the RDH or RDHF approaches on the extrapolation 
of the coupling constants at very low densities are observed |TB|, |T^. However, no detail 
procedures of the extrapolation were presented there . Since there is no direct restriction on 
the coupling constants at very low density , in Ref.|T^ the scalar and vector potentials of 
the DBHF results were extrapolated at low densities first by respecting their properties at 
nuclear matter. Then, the extrapolations of the coupling constants in the RDH or RDHF 
approaches are restricted by the scalar and vector potentials at the low densities. There- 
fore, the extrapolations of coupling constants at the low densities in different cases are on 
an equal level. As mentioned above, the contributions of isovector mesons are neglected in 
Refs. |T6|, p!7| . It is known that a realistic description of the nucleon-nucleon interaction in 
terms of meson exchange must include vr and p . Therefore , it is necessary to develop a 
relativistic theory for finite nuclei in the RHF with the isovector mesons vr , p included. Fritz, 
Miither [1^ discussed the vr contribution in the RDHF approach to the bulk properties of 
finite nuclei. They found that the inclusion of the vr-exchange terms in the RDHF slightly 
improve the agreement between calculation and experiment. On the other hand, Boersma 
and Malfliet PU| achieved a density dependent parametrization of the Dirac-Brueckner G 
matrix in nuclear matter which was called an effective DBHF . They used the effective DBHF 
to systematicly analyze a series of spherical nuclei. The results are in good agreement with 
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the experiments. It should be mentioned that the isovector vector meson p was not included 



in those calculations. In our previous brief report the effects of the isovector p meson 
in the RDHF on the bulk properties for finite nuclei are discussed but p tensor coupling was 
not included. It was found that the Fock exchange terms in the a-uj model reduce the charge 
radii, but have less influence on the binding energies. A large repulsion of vr contribution at 
the interior of nuclear is observed. As a result the energy levels of single particles becomes 
shallow at the presence of vr meson. Therefore the total binding energy is reduced and the 
charge radius is expanded. This effect is partly canceled when the p meson is included in the 
RDHF approach. But the contribution of tensor coupling of p meson has not been discussed 
in these works . In this paper, the effect of tensor coupling of p meson is included in the 
RDHF approach. The systematic study of finite nuclei in terms of the RDHF approach is 
investigated. 

The arrangement of this paper is as follows: The general formalism in this work is 
presented in Sec. 2 . The numerical results and main conclusions are included in Sec. 3 and 
Sec4. 

2 The formalism 

As in the one-boson-exchange(OBE) description of the NN interaction |p2||, our starting 



point is an effective Lagrangian density which couples a nucleon [ip) to two isoscalar mesons((j 
and a;) and two isovector ones( tt and p) with the following quantum number (J'^,T): 

a(0+,0),a;(l-,0),7r(0-,l),p(l-,l) 

and the electromagnetic field {A^) is also included. 

The effective Lagrangian density can be written as the sum of free and interaction parts: 

C = C, + Ci (1) 
The free Lagrangian density is given by 

>Co = ^{ilf^d^ - M)^ + ^{d^crd^a - ml) + ^m>^cu^ 
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with 



G^iy — dup^ — dfj^pu, 

where the meson fields are denoted hj a , , p^ and n , and nia , m^j , nip and are 
their masses , respectively. The nucleon field is denoted by ip which has a rest mass M . 
is the electromagnetic field. The interaction Lagrangian density is given by 

- — ^/'TsTm^^^ ■ (3) 

here r and are the usual isospin Pauli matrices. The effective strengths of couplings be- 
tween the mesons and nucleons are denoted by the coupling constants gi or fi {i = a,Lj, p, n) , 
respectively. Note that the pseudovector(PV) coupling for vrNN interaction is used , because 
the baryon self-energies are extremely large (about 40 times larger than their PV counter- 
part) at normal nuclear density if a pseudoscalar coupling is used , which has a drastic effect 
on the single-particle spectrum[^]. The present of tensor couplings makes the model La- 
grangian density no longer renormalizable and all physical observable should be calculated 
at the tree level 



2.1 Equations of motion 

The equation of motion for the meson fields are easily obtained from the Euler- Lagrange 
equation 



dC „ dC 

'd{d^(f)) 



— - = 0, (4) 
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with a meson field (j) . For instance , the a and u! fields are the solutions of 

{n + ml)(r = g^tljij, (5) 

(□ + ml)uj^ = gu^i^-f^i', (6) 

with the baryon current conservation {ip'') ^ip) = 0. Solving the equation for the meson 
fields , one then obtains 

(^ix) =g^ d^yD^{x - y)ip{y)ip{y), (7) 



uj'ix) =g^j d'yD>^^{x - y)i,{y)^,i,{y), (8) 

where -Do- (a; — y) and D^'^{x — y) is the a and uj meson propagator. Similar expressions are 
deduced for isovector mesons. 

Following standard techniques [p^], at the Hartree-Fock level , the average of Hamiltonian 



at the ground state can be written as 

< ^q\H\^o 51 / Ui{x)[-ia ■ V + 7om]f/„(x)c?x + / Ul{x)-ioT.H{x)Ua{x)dx 

- J Ul{x)-fo J i:Fix,y)Uaiy)dydx, (9) 
where Ua{x) is nucleon wave function and satisfies the orthogonality relation 

dyUi{y)Up{y) = Sap 

Only positive-energy states have been taken into account in the preceding derivation. 



2.2 Nuclear matter 

Because of the translational and rotational invariance in the rest frame of infinite nuclear 
matter and the assumed invariance under parity and time reversal, the nucleon self-energy 
produced by the meson exchanges in nuclear matter can , in general, be written as 

j:iK) = j:,iK) - ^oMK) + 1 ■ ks,(A;,) (lo) 
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where , Sq , Si, denote the scalar, time and space components of vector potentials, 
respectively. In general, they are functions of the four-momentum k^y of a nucleon and Fermi 
momentum kp. Based on the Feynman diagram rules one could derive the nucleon self- 
energy in nuclear matter. In the RHF approach, the isoscalar mesons in our Lagrangian 



density give rise to the following contributions to the self-energy [^, ^ 

= + r"rfggM(go)[^^e.(fc,g)-%50.(fc,g)], (11) 

7Tl/(j J-OTT rh Jo 

So(M = -i—rPB--^ rdqq[gle„{k,q) + 2glQ^{k,q)], (12) 



where 



{k + qY + 



Qi{k,q) = In 



{k — qY + ml 



_l_ _|_ lYi'^ 

" 4fcg ^®»(^' - 1' ^ = ^' ^' ^ 

k* = k(l + S,(A:,)), r = |k*| 

kl = ko + MK) = {k*^ + M*^f'\ 
The scalar and vector densities are 

2 i-f'F 



Ps = — I %^M{q)dq, 

TT^ Jo 



PB = ^,k% (14) 



2 

The contributions of the isovector mesons vr and p to the self-energy is given in the Appendix 
A. The first terms of the and Sq are the Hartree terms, which are energy- independent. 
The rest terms are the Fock terms, which are almost 1/k energy dependent. The coupled 
nonlinear integral equations (11-13) have to be solved self-consistently . 

With these nucleon self-energies , following the procedure discussed in Sec3.1, we can 
obtain the coupling constants and g^^ at each baryon density. 
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2.3 Finite nuclei 

At the case of spherical, closed-subshcll nuclei, a single-particle baryon state with energy 
is specified by the set of quantum numbers 

where q'q = — 1 for a neutron state and ^'q, = +1 for a proton state. The nucleon wave 
function can be written as 

Uo.{x)^-i f;l''\\y.Xi/2{q.). (15) 
r \ Fa{r)a -r J 

where Xi/2(?a) is an isospinor, and the angular and spin parts of the nucleon spinor can be 
written as 

ya{f) = J2 < Icy^'^aSaljarria > Yi'^" {r)xi/ 2 {s a) ■ 

{-I'd j^a 

The spinors Ua{r) are normalized according to 



The HF solution is obtained by requiring that the total binding energy 

E =< 0o|//|0o > -AM, 

is stationary with respect to variations of the spinors Ua{i-e- of Ga and Fa) such that the 
normalization relation is preserved 

S[E - Eaj Ui{r)Ua{r)d?r] = (16) 

a{occ) 

using the expression of nucleon spinor, after a lengthy derivation, the HF equations for the 
self-consistent wave functions [Ga, Fa) and energies Ea will be obtained. The radial Dirac 
equation take the following form: 

d_ f G„(r) \f - S^^,(r) M + Ea + Sf,„(r) - ^ 

dr [ F,{r) ) [M-Ea + Sf,„(r) + Eo^„ ^ + E|?_„(r) 
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here, E^^^, E^^ and are direct contributions to the self-energy and can be written as 

S?„(r) = [Ej(r) + E7W(r)]g„ (18) 

^U^) = (19) 
Kair) = E.(r) + [E^^(r) + E7(2)(r)]g. + 1(1 + g.)S,(r) (20) 
whereas and Ya come from exchange (Fock) contribution. The quantity Ka is 

{2ja + l){la-ja)- 

In this work we consider nuclei with a closed proton and neutron shell only, therefore, 
the isovector pseudoscalar meson yield no contributions in the Hartree approximation. The 
Hartree contributions come from a, cu, p mesons and Coulomb force are given as 

Ea(r) = -ga{pB{r))m^ / g^{pB{r'))ps{r')Io{mar<)KQ{m„r>)r''^dr' (21) 

Ew(r) = 9Lu{pB{r))maj / gw{PB{r'))pB{r')Io{mu,r<:)Ko{m^r>y^dr' (22) 

Jo 

E^(r) = gimp / [pbA^') - PB,nW)M^P^<)Mmpr^y^dr' (23) 

t/ 

= -f^™. /„°°l^'-.>'(''') - PT,„(r')l[^/„(m,r<)A'„(m,r>)]r'2<ir' (25) 

Sj(r) = -i^fH I [pT,p(r') - pT,n(r')]/i(m,r<)iri(m,r>)r'^dr' 

-[pT,p(r)-pT,n(r)]} (26) 

with the definitions 

P5„=^ E jt[G^(r)-F,^(r)] 

b{norp) 



pT,norp = 3fo[2G'fe(r)Ffc(r)] 

b{norp) 

and r< (r>) is the smaller (larger) of the r' and r. 

The functions and i^^L(3:^) arise from the multiple expansion of the meson propagator 
in the coordinate space and are defined using the modified spherical Bessel functions of the 
first and third kind I and K: 

The explicit expressions of exchange contributions Xa and Yq, are given in the Appendix B. 



3 Results and Discussions 
3.1 Parameterization of DBHF 

The self-energy obtained in the RHF is of very weak energy-dependence and too strong 
density-dependence due to the fact that the short range correlation has not been considered 
. Attempts to incorporate the effects of the short range correlation described in the DBHF 
approach have been made by introducing density- and momentum-dependent effective cou- 
pling constants of mesons in the RHF approach . In order to make comparison with the 
DBHF results, the scalar and vector potentials can be obtained by 

Us{k, kp) = — , Uo{k, kp) = — ^-p^ , (27) 

where — \Jk*'^ -\- M*^ — Eo- The momentum dependence of the potentials usually are 
relatively weak and neglected in the description of ground state properties of finite nuclei . 
Therefore , the momentum average within the Fermi sea is performed, 

At very low density of nuclear matter the DBHF results are not reliable and remain un- 
known. Therefore, the extrapolation of the coupling constants outside the density points in 
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the DBHF has to be done when they are apphed to the calculation of finite nuclei. In order 
to remove the sensitivity, the extrapolation of scalar and vector potentials and Uo of the 
DBHF results at low densities are done by setting = 0, Uq = at p = 0. It is known 
that the scalar and vector potentials in RMF or RHF are almost linear dependent on the 
density. Due to the two-body correlation the scalar and vector potentials approach to zero 
smoothly as the density goes to zero. A polynomial fit of the scalar and vector potentials 
with respect to the density are performed. The extrapolation and interpolation of and 
Uo are shown in Fig.l, where the circles are the DBHF results in nuclear matter using the 
Bonn A potential [|13|. The density dependence of the coupling constants are then adjusted 
in the cases of RMF or RHF with or without isovector mesons to reproduce the nucleon 
self-energies at each densities resulting from the DBHF. The nucleon and a and uj meson 
masses are chosen to be the same as the DBHF calculation, where M = 938.9 MeV, m^. = 
550 MeV, m^j = 782.6 MeV. The pseudo-vector coupling for vrNN and vector and tensor 
coupling for pNN are adopted. The masses and coupling constants of isovector mesons are 
fixed to be m^ = 138 MeV, m^ = 770 MeV, g = 0.08, || = 0.55 , ^ = 3.7 §. The density 
dependence of coupling constants in different cases are shown in Fig.2. The presence of the 
pion introduces a large repulsive force, so the scalar coupling constant becomes larger and the 
vector coupling constant gets smaller to balance the repulsive force, especially at normal and 
high densities. However, the pion contribution is partly canceled by the presence of vector 
part of p meson at the symmetric nuclear matter . The tensor coupling of p meson has large 
effect at high density . As a result , the coupling constant of a meson becomes larger than 
that of uo meson at high densities . The results obtained in this paper are somewhat different 



from those in ref. The reason for this discrepancy is that the zero-range components of 
the pion-exchange were removed there. The cases with and without the contact interactions 
have been discussed in more detail in ref. [^]. It is found in our calculations that the effects 
of the contact interactions mainly cause a renormalization in g^j and coupling constants. 
The removal of the zero-range components of the pion and rho exchange would increase the 
binding energy and reduce the charge r.m.s radius. No qualitative improvement has been 
found. Therefore, only the cases with the zero-range components of the vr and p exchange 
are presented in this paper. 
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3.2 The ground state properties of finite nuclei 

The ground state properties of four stable doubly closed-shell nuclei ^^O, ^''Ca, ^^Ca and 
^°Zr are calculated with these density dependent coupling constants in the RDH and RDHF 
approaches . The set of coupled differential equations(17) is solved in the coordinate space 
following the method of Ref. P]. The self-consistence is achieved by an iterative procedure. 



It is different from a matrix diagonalization method adopted in Refs. |1T6|, p!7| , where a ter- 
mination of a complete set of bases has been performed both for baryon and meson results. 
Our computer code has been carefully checked with the results of ref.0, ^ . 

In order to investigate the effect of the density-dependence, the results for ^^O and ^°Ca 
obtained in the RMF and RHF with of a + (RHFl) , a + u; + vr (RHF2) and a + uj+ 
77+ p(RIIF3) are list in Table 1. The coupling constants are determinated to reproduce the 
DBHF results in nuclear matter (OBE potential A) at saturation density Kp = 1.40 fm~^. 
It should be mentioned that the results are different from those of the usual RMF and RHF 
calculations , where the coupling constants and scalar meson mass are adjusted to reproduce 
the empirical saturation properties of nuclear matter as well as the rms charge radius of ^°Ca 
. Because of relative large saturation density obtained in the DBHF and the scalar meson 
mass rricr = 550 Mev adopted in this calculation , the binding energies and rms charge radii 
calculated here are both much smaller than the experimental data . However , the main 
purpose of Table 1 is to show the difference in various cases mentioned above as well as the 
effect of the density-dependence in comparison with Table 2. The calculations with density- 
dependent effective interaction are performed, where the coupling constants at each baryon 
density come from the parametrization of the DBHF result in nuclear matter as discussed 
in section 3.1. Various cases, RDH, RDHF with cr + cj(RDHFl), a + cu +7i (RDHF2) and 
cr +u! + TT + p (RDHF3) are investigated and the results for the nuclei ^^0,^°Ca , ''^Ca 
and ^°Zr are displayed in Table 2 . In order to investigate the sensitivity of the p meson 
couphng constant , two values of p coupling constants are adopted in the calculations : ^ = 
0.55(RDHF3A) and 0.99(RDHF3B) without tensor coupling. The results with pNN tensor 
coupling (|| = 0.55, ^ = 3.7) are given as RDHF3C. 

The importance of the density dependent approaches is clearly demonstrated in Table 
1 and Table 2. The calculations in either RMF or RHF with constant coupling constants 
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RMF 


RHFl 


RHF2 


RHF3 


160 










EB/A{MeY) 


-5.62 


-6.02 


-4.86 


-5.67 


rc/i(fm) 


2.48 


2.39 


2.53 


2.57 


lsi/2(Mev) 


44.34 


45.08 


40.78 


43.21 


lp3/2(Mev) 


18.96 


21.15 


17.50 


18.51 


l]9i/2(Mev) 


9.62 


7.85 


9.49 


10.98 












EB/AiMeY) 


-6.36 


-6.69 


-5.84 


-6.38 




3.14 


3.04 


3.16 


3.22 


14/2 (Mcv) 


16.54 


18.61 


15.72 


16.35 


2si/2(Mev) 


7.07 


5.60 


7.98 


8.81 


14/2 (Mev) 


6.92 


5.70 


6.50 


7.79 



Table 1: Ground state properties of ^^O and ^°Ca calculated by RHA and RHF. The binding 
energy per nucleon Eb/A, the charge rms radius Tc and single-particle energies of proton 
states. 
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RDH 


RDHFl 


RDHF2 


RDHF3A 


RDHF3B 


RDHF3C 


Exp. 


160 

EB/A{Mew) 
rc/i(fm) 


-7.44 
2.59 


-7.48 
2.50 


-6.96 
2.64 


-7.29 
2.61 


-7.57 
2.59 


-7.41 
2.68 


-7.98 
2.73 


lsi/2(Mev) 
lp3/2(Mev) 
lpi/2(Mev) 


43.97 
21.77 
16.16 


43.87 
23.60 
16.08 


41.11 
21.11 
15.60 


42.76 
21.85 
16.07 


44.13 
22.47 
16.47 


42.98 
21.31 
15.72 


40 ± 8 
18.4 
12.1 


EB/A{Mev) 
rcft(fm) 


-7.88 
3.26 


-7.91 
3.17 


-7.49 
3.29 


-7.74 
3.27 


-7.94 
3.26 


-7.81 
3.35 


-8.55 
3.48 


14/2 (Mev) 
2si/2(Mev) 
lo?3/2(Mev) 


19.27 
13.69 
13.29 


21.17 
14.08 
13.48 


19.08 
14.20 
12.77 


19.52 
14.19 
13.09 


19.89 
14.14 
13.35 


18.95 
13.78 
12.67 


15.5 
10.9 
8.3 


EB/A(Mev) 
rcfe(fm) 


-8.02 
3.27 


-7.96 
3.17 


-7.45 
3.30 


-7.59 
3.28 


-7.70 
3.27 


-7.60 
3.37 


-8.67 
3.47 


14/2 (Mev) 
2si/2(Mev) 
14/2 (Mev) 


24.35 
17.25 
18.68 


29.17 
19.77 
21.79 


24.34 
19.79 
23.05 


26.37 
20.90 
24.92 


27.96 
21.77 
26.38 


26.72 
20.61 
23.53 


20.0 
15.8 
15.3 


90 Zr 

EB/A{MeY) 


-7.94 
4.00 


-7.92 
3.89 


-7.55 
4.02 


-7.67 
4.00 


-7.78 
3.99 


-7.68 
4.10 


-8.71 
4.27 


2p3/2(Mev) 
l/5/2(Mev) 
2pi/2(Mev) 


11.01 
13.14 
9.42 


13.09 
15.61 
11.16 


12.86 
16.67 
11.77 


13.61 
18.12 
12.49 


14.19 
19.24 
13.05 


13.52 
16.93 
11.91 


11.0 
12.3 
9.5 



Table 2: Ground state properties of ^^0 , "^^Ca , ""^Ca and '*Zr calculated by RDH and 
RDHF. The binding energy per nucleon the charge rms radius r^ and single-particle 

energies of proton states. 
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produce much smaller binding energies of nucleon in comparison with the experiments. In 
contrast, the density-dependent interactions increase both binding energy and charge radius, 
which imply the removal from the so-called Coester bandQ. The results in the relativistic 
density-dependent approaches are largely improved and closer to the experimental values. 



The slight differences from ref. [|1^, are due to the different extrapolation procedures. 
The Fork exchange term in the a-uj model reduces the charge radii, but has less influence 
on the binding energies. A large repulsion of pion contribution at the interior of nucleus 
is found. As a result, the energy levels of single particles become shallow at the presence 
of pion. Therefore, the total binding energy is reduced and the charge radius is expanded. 
This effect is partly canceled by the p meson exchange contribution. The charge radii of 
nuclei calculated in the RDHF with isovector mesons are much close to those obtained in 
the RDH, which can also be observed in the charge density distributions. In comparison of 
the RDHF3A and RDHF3B , the results are not sensitive to the strength of the p meson 



,2 



coupling. The binding energies for the strong coupling constant of p meson ^ = 0.99 is 
about 2% bigger than those for ^ = 0.55 and the rms radius is reduced less than 1%. With 
p tensor coupling ,it can be found that the results of both binding energy and rms charge 
radius are improved . The binding energies of nuclei is increased slightly, but the charge 
radii of nuclei are improved largely in comparison with experiments. 

Figure 3 shows the density distribution of nuclei. The dash-dotted curves are the results 
of the RHF3, corresponding to the fourth column in Table 1. The dashed, dotted and solid 
ones correspond to those of the RDH and RDHF without and with isovector mesons (RDHFl, 
RDHF3A, RDHF3C) , respectively . It is found that the charge densities are reduced at 
the nuclear interior and have long tail due to the relative strong coupling constants at the 
nuclear surface in the density dependent calculations. The Fock contribution of a and uj 
in the RDHF produce a squeezing effect and give a large central density. The repulsive 
contribution of vr reduces the interior density and the results of the RDHF3A are very close 
to those of the RDH . Though the density at center is still higher than the experimental 
results, it is a parameter-free calculation in the sense that no parameters are adjusted for 
the calculations of the many-body problem. However, the results with density dependent 
coupling constants are in a reasonable good agreement with the experiments. 

As well known , nonrelativistic BHF calculations with various tow-body nucleon-nucleon 
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potentials, such as Reid soft-core ,Hamada- Johnson potentials , reveal a "Coester" band in 
dependence of 1/rc on Eb/A. In the RMF or RHF calculations, the dependence of l/vc on 
Eb/A at a fix Kp, with variation of the scalar meson mass and therefore the variation of the 
coupling constants , formed a new "Coester" band. A better estimate of the merits of the 
present work can be expected upon the comparison with the BHF " Coester" band and RHF 
"Coester" band. Those "Coester" bands are plot in Fig.4 for ^^O and ^°Ca , the dash-dotted 
line indicates the BHF "Coester" band, which is taken from a so called generalized BHF 
calculations by Kiimmel et al.[^. The dotted and dashed lines represent the RHF results , 
which are obtained in the RHF3 by varying the scalar meson mass as well as the different 
coupling constants to reproduce the nuclear matter saturation properties at Kp = 1.40fm~^ 
(1) and Kp = 1.30fm~^ (2) resulting from the DBHF approach , respectively . The results 
in the case of the RDHF3C are displayed by a solid line in the figure. The density dependent 
approach forms a new line away from all of the "Coester" band of conventional BHF and 
seems to be much closer to the experimental values than the BHF and RHF calculations. 

The spin-orbit splittings of nuclei are given in Table 3. It can be seen that the spin- 
orbit splitting in the RMF and RHF approaches is larger than the experimental data, which 
indicates the larger spin-orbit force. It is known that the spin-orbit force is related to the 
derivative of the potentials with respect to the space and is a surface effect. The density 
dependent approaches reduce the sharp surface and, therefore, reduce the spin-orbit splitting. 
The Fock terms of the a and uj exchange increase the spin-orbit splitting, while the vr and p 
exchanges give the opposite contribution. A large reduction of the spin-orbit splitting due 
to pion-exchange is found, especially for heavy nuclei with large neutron excess. A flip of the 
spin-orbit splitting in ^osp]^ jg observed in the calculation of RDHF with all mesons included, 
which is certainly not physical. It might indicate that the free coupling constants of the 
isovector mesons adopted in the RDHF are too strong for the nuclear structure calculation. 
A similar observation was obtained in the calculation of the relativistic optical potential in 
the RHF approach]^. The density dependent coupling constants for isovector mesons may 
also be required. 

It is known that the single particle densities are not directly provided by experiments . 
The only way to gain some insight in the single particle distribution is to study the difference 
between density distribution of nearby nuclei. In Fig. 5 ,we give the charge distribution 
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Nuclei RMF RHFl RHF3 RDH RDHFl RDHF3A RDHF3C Exp. 



160 


9.34 


13.30 


7.53 


5.61 


7.52 


5.78 


5.60 


6.3 




9.62 


12.82 


8.56 


5.78 


7.69 


6.43 


6.28 


7.2 


48Ca 


9.32 


12.52 


5.32 


5.67 


7.37 


1.45 


3.19 


4.3 



Table 3: Spin-orbit splittings of protons for the Ip shell in ^^O and the Id shell in ^°Ca and 
^^Ca. 



Nuclei 


RHF3 


DBHF 


RDHFl 


RDHF3C 


Exp. 


160 


-0.03 


-0.03 


-0.03 


-0.02 


-0.02 


^"Ca 


-0.05 


-0.06 


-0.05 


-0.04 


-0.07- 0.10 


^^Ca 


0.23 


0.13 


0.15 


0.22 


0.16 - 0.23 


902r 


0.11 


0.04 


0.07 


0.11 


0.07 



Table 4: Neutron skin thickness A„p = r„ - rp for ^^O , ^°Ca , ^^Ca , ^^Zr . The DBHF 
results are taken from [20]. 

difference between ^°Ca and ^^Ca, multiplied by . The difference of the neutron densities 
between ^°Ca and ^'^Ca are plotted in Fig. 6 . The shaded area presents the experimental 
data, the dashed curve is the results of the RDH , the dotted one the results of RDHFl , and 
the solid one is obtained in the case of RDHF3C. It can be seen that the results obtained in 
the RDHF3C are superior to those of the RDH and RDHFl in comparison with experiments. 
It means that the isospin dependence can not be correctly described by the RDH as well as 
RDHFl without including isovector mesons. 

Neutron skin thickness is an important quantity to study isotope shifts. It is defined 
as the difference between neutron and proton rms radii: A„p = r„ — Tp. The neutron 
skin thickness of ^^O , ^°Ca , ^^Ca and ^°Zr are given in Table 4 , and the A„p versus the 
asymmetry parameter (N-Z)/A for ^''Ca, ^'^Ca and ^°Zr are shown in Fig. 7 . The results 
of the RDHF3C seem to be similar to those of the RMF , which are different from what 
obtained in ref. PU| . More information of of isospin dependence the ground state properties 
are required . 
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4 Conclusion 



In summary , the RDH and RDHF approaches with the density-dependent effective coupling 
constants of isoscalar mesons can incorporate the DBHF results and contain the nucleon- 
nucleon correlation effects . Inclusion of the NN correlation let to a substantical improvement 
in the microscopic description of bulk properties of nuclei . The Fork exchange terms are not 
neghgible , though the exchange contributions are relatively weak than those of the Hartree 
direct term in the relativistic approach and their contribution to the binding energy may be 
compensated by the variation of the coupling constants . The important contributions from 
the isovector meson tt and as well as , to some extent , p meson are not included in the mean 
field approach . It is found that the isovector mesons tt and p play an important role in the 
spin-orbit splitting as well as the isospin dependent quantities . The tensor coupling of p 
meson gives a constructive contribution to the bind energy , especially the rms charge radius 
, therefore improves the agreement with the experimental data . More information of the 
isospin denpendence of nuclear properies is required to provide constraints on the coupling 
constants of isovector mesons tt and p in the nuclear medium . 
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Appendix A. Self-Energy in Nuclear Matter 



Base on Hartree-Fork approach , the nucleon self-energy in nuclear matter coming from 
the contributions of isovector pseudoscalar meson vr can be written as follow: 

= J,'- dqq{J^fM{2kq - X©.) (^-1) 

= J/;. dqq{^n-lmlQ^ + 2kq] {A.2) 

K(k) = /o^- dqqi^nqQke^ - Q{e + q')^^] {A.3) 

The contribution coming from isovector vector meson p can be write as : 



= /o " dqq{-2glMe, + 3{I^fM{2kq - Im^G, 
+3^,(4)[gge,-2A;Q$,]} 



{AA) 



^oik) = /o " dqq{gle, + {^n2kq - im^Qj} (^.5) 
= lo' dqq{2glQ^, + 2{^f[kqQQ, - Q{k^ + _ 1^2)$^ 



~3g,{^){kMep-2qM^p)} 
Appendix B. Fock Term Expressions 



(A.Q) 



The quantities X and Y of eq.|I^ can be written as the sum of contributions coming 
from different mesons . In the following , we often need the reduced matrix elements of the 
tensorial operators Y[^{r) and 



mk 



They are given as follow: 



,„ , (47r)-ijJ,L(-l)-'-^-^-H ^\ ^ I iila + h + Liseven 
< a Yr 6 >= <^ I 1 _1 ' 



2 2 



< a\\TjL\\b >= mH-lY'^jajbUjL 



la L If, 
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where J = V2J + 1- 

The contribution coming from scalar moson a is determined to be 

;:;:f') - ..(PB«)".j-i:*..(2:',)i:i<«iinii^>r 

X / g,{pB{r'))\G,G^-F,F^\,U{m,r^)Ki^{m,r^)dr' , (5.2) 

the sum over b running over occupied states . 

The expressions for the vector meson uj are sphted into timehke and spacehke parts , due 
to the respective 70 and 7 couphng. The time component is 

X / 5c.(PB(r'))[GaG'„ + F,Fb],,/i(m,r<)i^z,(m,r>)rfr' , (S.3) 

JO 

The space component is 



M (..4) 



2 



2 



I, ^ I, J \l -I Q 

X lL{m^r<)KL{m^r>)dr' 

where a' = (g^, n^, j^) with l'^ = 2ja - k ■ 

The contribution comes from isovector pseudoscalar meson tt with pseudovector couphng 
is as follows 



/ [{i^ab + a{L2)}GaGb - {Kab - Q!(L2)}F„Fb]^/i?L ^^(m^r, my)dr'} 
Jo 



L LiM 

roo 

X 
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where we introduced the notation 

f -L iiL, = L-l, 
a{Li) = { 

[ L + 1 if Li = L + 1, 

-RliL2(^^) ^^') — lLi{'m'r')KL^{mr')0{r' — r) + KL^{mr)lL2{jnf'')9{r — r'). 

The Li and L2 can only take two values L + 1 or L — 1 . 

The vector part of the pNN coupling in our Lagrangian gives X{r) and Y{r) which are 
formally idential to those of the uj meson , expect for the isospin factor 8q^q^ replaced by 2 
- 5q^q^ and , gu}{pB{f)) replaced by nip , Qp . The tensor term gives rise to two types of 
contributions . They are proportional to and fpQp , and they are denoted , respectively 
, by {X^'^\Y^'^^) and F^^^^) . Furthermore , they can be splited into timehke and 

spacelike component . 

The time component of X^"^^ and Y^'^^ is written as 



rup 



L LiM \ Fb(r)[Kab + Ol(Li)] J 

POO 

X / [{iiab + Oi{L2)}GaFh- {kab- a{L2)}FaGb]r'RLi,L2i^P^^^P^')'^^'} ■ 
J 

where kab ^ Ka - Kb . 

The space component of X^^^ and Y^"^^ is 

(< a'WTjr \\b' > \ 
/ [<a\\TjL2\\b>GaGb-<a'\\TjL2\\b'>FaFb]r' 
<a\\TjL,\\b> J Jo 

[mpRLiL2 (rupr, nipr') - ^^^^^y^]dr', {B.l) 



where we have introduced 

hj{Li)=LL, 





L 






L 













1 


^1 
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The time components of VT contributions are 



-Gbir)<a'\\TLL,\\b> 

F, {r) < a\\TLLA\b' > 

/ [< a\\YL\\b > GaGb+ < a\\YL\\b' > FaF,]r>SLL,{r,r')dr' (5.8) 
Jo 

-Fb{r) < a'\\YL\\b' > 

G, (r)<a\\YL\\b> 

POD 

/ [< a\\nL,\\b' > GaF,+ < a'\\TLL,\\b > F„G',],,5LL,(r', r)dr'}. 
Jo 

The space component of VT contribution are 



{ 



X 



+ 



X 



{ 



Fb{r) < a'\\TjL,\\b' > 
Gb{r) < a\\TjL,\\b > 



roo 

X / [< a\\TjL\\b' > GaFk- < a'\\TjL\\b > FaGh]r'SLLAr:r')dr' (5.9) 
Jo 



+ 



X 



Gt,{r)<a'\\TjL\\b> 
F,{r) < a\\TjL\\b' > 



roo 

/ [< a\\TjL,\\b > GaGb- < a'\\TjL,\\b' > FaF,]r'SLL,{r,r')dr'}. 
Jo 



where 



SLLi{r,r') = lL^{mpr)KL{mpr')9{mp{r' - r)) - lL{mpr')KL^{mpr)9{mp{r - r')). 
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Figure captions 



Figure 1. Scalar and vector potentials and Uq as functions of the density in nuclear 
matter. The circles are the DBHF results using Bonn A potential[l]. The curves are 
obtained in terms of interpolations and extrapolations. 

Figure 2. Density dependent coupling constants of cr and uj, g„ and g^. They are deduced 
by reproducing the scalar and vector potentials of the DBHF results at each density 
from RMF analyses (a) and RHF analyses (b). The solid and dotted curves in (b) are 
corresponding to the cases of 0"+ uo and a+cu+vr+p in RHF. The circles on curves refer 
to the density points of the DBHF results. 

Figure 3. Charge density distribution of various nuclei . The cureves are the results of 
RDH (dashed one), RDHF with a + uj only (dash-dotted) and RDHF with a + uo + 
IT + p { solid for 3C and dotted for 3 A ) and RHF with a + u + tt + p (dence-dotted 
one ) . 

Figure 4. Binding energy versus l/vch for ^^O and ^°Ca . The dash-dotted lines is taken 
from the work of Kiimmel et al. The dotted and dashed lines are obtained by 
RHF3 with changing a meson mass for Kp = 1.4 fm^^ (1) and Kp = 1.3fm~^ (2) , 
respectively . The solid lines are the results of RDHF3C . The experimental data are 
displayed by a star . 

Figure 5. Difference between charge densities of ^"Ca and '^^Ca multiplied by . The 
dashed curve corresponds to the RDH , the dotted curve to the RDHFl . The solid 
line corresponds to the RDHF3C . The shaded area indicates the experimental data. 

Figure 6. Same as Fig. 4 for difference between neutron densities of ^'^Ca and ^°Ca. 

Figure 7. A„p versus the asymmetry parameter for various cases and the DBHF approach 
taken from ref. |20[ 
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